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The harmonic oscillator with dissipation is studied within the framework of the Lindblad 
Q\ • theory for open quantum systems. By using the Wang-Uhlenbeck method, the Fokker-Planck 

equation, obtained from the master equation for the density operator, is solved for the Wigner 
| distribution function, subject to either the Gaussian type or the 5-function type of initial 

conditions. The obtained Wigner functions are two-dimensional Gaussians with different 
■ widths. Then a closed expression for the density operator is extracted. The entropy of 

the system is subsequently calculated and its temporal behaviour shows that this quantity 
^ ; relaxes to its equilibrium value. 

D 

& ■ 1 Introduction 

> '■ 

In the last two decades, the problem of dissipation in quantum mechanics, i.e. the consistent 
description of open quantum systems, was investigated by various authors p], |], ||. Because 
dissipative processes imply irreversibility and, therefore, a preferred direction in time, it is gen- 
erally thought that quantum dynamical semigroups are the basic tools to introduce dissipation 
in quantum mechanics. In the Markov approximation the most general form of the generators 
of such semigroups was given by Lindblad 0. This formalism has been studied for the case of 
damped harmonic oscillators |7| || || and applied to various physical phenomena, for instance, the 
damping of collective modes in deep inelastic collisions in nuclear physics |TUJ and the interaction 
of a two-level atom with the electromagnetic field |TT] . 



In the present work, also dealing with the damping of the harmonic oscillator within the Lind- 
blad theory for open quantum systems, we will explore the physical aspects of the Fokker-Planck 
equation which is the c-number equivalent equation to the master equation for the density opera- 
tor. Generally the master equation gains considerably in clarity if it is represented in terms of the 
Wigner distribution function which satisfies the Fokker-Planck equation. It is worth mentioning 
that these master and Fokker-Planck equations agree in form with the corresponding equations 



formulated in quantum optics [12 , 13 , \L4\ , |15[ |16| . 



The content of the paper is arranged as follows. In Sec. 2 we review the derivation of the 
master equation of the harmonic oscillator. In Sec. 3 we transform the master equation into the 

Then the Fokker-Planck 
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Fokker-Planck equation by means of the well-known methods [|17], |TI 
equation for the Wigner distribution, subject to either the Gaussian type or the 5-function type 
of initial conditions, is solved by the Wang-Uhlenbeck method. Sec. 4 derives an explicit form of 
the density operator involved in the Lindblad master equation, formulates the entropy using the 
explicit form of the density operator and discusses its temporal behaviour. Finally, concluding 
remarks are given in Sec. 5. 



2 Master equation for the damped harmonic oscillator 

The rigorous formulation for introducing the dissipation into a quantum mechanical system is that 
of quantum dynamical semigroups 0, |3], According to the axiomatic theory of Lindblad |J, 
the usual von Neumann-Liouville equation ruling the time evolution of closed quantum systems 
is replaced in the case of open systems by the following equation for the density operator p: 

d ^ = L(Up))- (1) 

Here, $t denotes the dynamical semigroup describing the irreversible time evolution of the open 
system in the Schrodinger representation and L the infinitesimal generator of the dynamical semi- 
group Using the structural theorem of Lindblad [0] which gives the most general form of the 
bounded, completely dissipative Liouville operator L, we obtain the explicit form of the most 
general time-homogeneous quantum mechanical Markovian master equation: 

^ = L( P (t)) = -t[H,p(t)} +^£([^(t),^ + ] + [Vj,p(t)V+]). (2) 
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Here H is the Hamiltonian of the system and the operators Vj and are bounded operators on 
the Hilbert space of the Hamiltonian. 

We should like to mention that the Markovian master equations found in the literature are 
of this form after some rearrangement of terms, even for unbounded Liouville operators. In this 
connection we assume that the general form of the master equation given by (2) is also valid for 
unbounded Liouville operators. 

In this paper we impose a simple condition to the operators H, Vj, that they are functions 
of the basic observables q and p of the one- dimensional quantum mechanical system (with [q, p] = 
ih) of such kind that the obtained model is exactly solvable. A precise version for this last 
condition is that linear spaces spanned by first degree (respectively second degree) noncommutative 
polynomials in q and p are invariant to the action of the completely dissipative mapping L. This 
condition implies J5J that Vj are at most first degree polynomials in q and p and H is at most a 
second degree polynomial in q and p. Then the harmonic oscillator Hamiltonian H is chosen of 
the form 

H = H + ^(qp + pq), H = —p 2 + —q 2 . (3) 



With these choices the Markovian master equation can be written ||: 

dp % % i 

-j t = - fr[H ,p\ - t^(X + p)[q, pp + pp] + ^(A - p)\p, pq + qp] 

Dpp -[q, [q,p]} - ^[p, \p,p}} + £g([q, [p, p}} + [p, [q,p]}), (4) 



where D pp , D qq and D pq are the diffusion coefficients and A the friction constant. They satisfy the 
following fundamental constraints ||: 

i)D pp >0, ii)D qq >0, iii) D pp D qq — D 2 q > - . (5) 
In the particular case when the asymptotic state is a Gibbs state 

pc(oo) = e^kr /Tie~kr , (6) 

these coefficients reduce to 

„ X + p, L , tiuj A — p h J _ 1 hw n 

i- ) rm = ftmu; coth — — -, JJ ao = coth — — , D vo = 0, (7) 

where T is the temperature of the thermal bath. 

3 Wigner distribution function 

One useful way to study the consequences of the master equation (4) for the density operator 
of the one- dimensional damped harmonic oscillator is to transform it into more familiar forms, 
such as the equations for the c-number quasiprobability distributions Glauber P, antinormal 
ordering Q and Wigner W associated with the density operator [30]. In this case the resulting 



differential equations of the Fokker-Planck type for the distribution functions can be solved by 
standard methods WA [Bl ETI] employed in quantum optics and observables directly calculated as 



correlations of these distribution functions. 

The Fokker-Planck equation, obtained from the master equation and satisfied by the Wigner 
distribution function W(x\,X2,t) of real variables X\,x<i corresponding to the operators q,p 



has the form [20[ 



where 



muj 1 . . 



I _f' A -M -w \ n w _ 1 ( mujD qq D pq . 
" l u X + p)' ^ ~ h\ D pq Djmw)- 



Since the drift coefficients are linear in the variables X\ and x 2 and the diffusion coefficients are 
constant with respect to x\ and x 2 , 



Eq. 



Following the method developed by Wang and Uhlenbeck [£3 



(9) describes an Ornstein- Uhlenbeck process [^2, p3|] . 

we shall solve this Fokker-Planck 



equation, subject to either the wave-packet type or the 5-function type of initial conditions. 

f ) When the Fokker-Planck equation is subject to a Gaussian (wave-packet) type of the initial 
condition (xw and x 2 o are the initial values of x\ and x 2 at t = 0, respectively) 

,(xx, x 2 , 0) = — exp{-2[(x! - x w ) 2 + (x 2 - x 20 ) 2 ]}, (11) 



irh 



the solution is found to be 



W w (x 1 ,x 2 ,t) = —. exp{--^-[0 w (xi - xi) 2 + 4> w (x 2 - x 2 ) 2 + Xw(xi - X\){x 2 - x 2 )]}, (12) 

nhu!\ \B W \ 



where 

B w = 9l g 2 - \gl 9l = Ql = ^e 2Ai + ^(e 2At - 1), g 3 = 2[e~™ + ^(1 - (13) 

<t> w = gia* 2 + 92a 2 - 93, ipw = 9i + 92 ~ 93, Xw = 2(g x a* + g 2 a) - g 3 (a + a*). (14) 

We have put a = (/i — iQ)/u,A = —X — ifl and d\ = (a 2 mu)D qq + 2aD pq + D pp /mu)/h, d 2 = 
[mujDqq + 2fiD pq /uj + D pp /muj)/Ti and fi 2 = uj 2 — fi 2 . The functions x± and x 2 , which are also 
oscillating functions, are given by 

Xi = e~ xt [xio(cos fit + ^ sin fit) + ^20^ sin fit], (15) 

x 2 = e _A * [220 (cos fit — — sin fit) — ^10^ sin fit]. (16) 

2) If the Fokker-Planck equation (9) is subject to the 5-function type of initial condition, the 
Wigner distribution function is given by 

W(x u x 2 ,t) = ^^exp{-^[0 d (xi - xi) 2 + ip d (x 2 -x 2 ) 2 + Xd(xi ~ Xi)(x 2 -x 2 )]}, (17) 

wfiujx \B\ B 



where 



B = hf 2 -fl fi = f; = ^(e 2At -l), f3 = j(l-e- 2Xt ) 



<Pd = fia* 2 + f 2 a 2 -2f 3 , Vrf = /i + /2-2/ 3 , Xd = 2[/ia* + f 2 a - / 3 (a + a*)]. (19) 

So, the Wigner functions are 2- dimensional Gaussian distributions with the average values X\ and 
x 2 and different widths. 

When time t — > 00, x% and x 2 vanish and we obtain the steady state solution: 

W(x u x 2 ) = 1 exp[-~ £ )■>:■>■/.. (20) 



27r v /detcr iy (oo) 2 ij=i, 2 

The stationary covariance matrix a w (oo) can be determined from the algebraic equation 

Aa w (00) +/ (00) A T = Q W . (21) 



4 Entropy and effective temperature 



Entropy is a quantity which may be visualized physically as a measure of the lack of knowledge 
of the system. When we denote by pit) the density operator in the Schrodinger picture for the 
harmonic oscillator, the entropy S(t) is given by 

S(t) = -kTr(p\np). (22) 

For calculating the entropy we shall compute straightway the expectation value of the logarithmic 
operator < lnp >= Tr(plnp). Accordingly, the problem amounts to derive the explicit form of 
the density operator for the damped harmonic oscillator. 

To get the explicit expression for the density operator, we use the relation p = 2TrhN{W s (q,p)}, 
where W s is the Wigner distribution function in the form of standard rule of association and N 



is the normal ordering operator |17|, |24] which acting on the function W s (q,p) moves all p to 
the right of the q. By the standard rule of association is meant the correspondence p m q n — > 
q n p m between functions of two classical variables (q, p) and functions of two quantum mechanical 
canonical operators (q,p). The calculation of the density operator is then reduced to a problem of 
transformation of the Wigner distribution function by the N operator, provided that W s is known. 
A special care is necessary for the N operation when the Wigner function is in the exponential form 
of a second order polynomial of q and p. The Wigner distribution function previously obtained 



corresponds however to the form of the Weyl rule of association [p5|| . The solution (12) of the 
Fokker-Planck equation (9), subject to the wave-packet type of initial condition (11) can be written 
in terms of the coordinate and momentum as: 

W(q,p,t) = ——^exp{--[(j)(q- < q >) 2 + ^(p- < P >? ~ 2x(q- < q >)(p- < V >)]}, (23) 
where 



2h 



< q >= \ Xx, < p >= v2hmujx 2 , (24) 



hu 2 1 
4D? muj 



= CT PP =< q 2 > - < q > 2 = VV,, , (25) 



Tiuj 2 

^ = Oqq =< P 2 > ~ < V > 2 = -^2 mUJ 4>w, (26) 

X = <Tpq(t) = ~ < qp + pq > - < q >< p >= ^fiXw, 5 = - \ 2 (27) 

and < A >= Tr(pA) denotes the expectation value of an operator A. The Wigner distribution 
function (23) can be transformed into the form of standard rule of association [2£| by 



1 d 2 

W s (q,p) =exp(-ih-^-g-)W(q,p). (21 



Upon performing the operation on the right-hand side, we get the Wigner distribution function 
W s , which has the same form as the original W multiplied by h but with \ ~ ifr/2 in place of x- 
The normal ordering operation of the Wigner function W s in Gaussian form can be carried out by 
applying McCoy theorem |24], |27], ESI . The explicit form of the density operator is the following: 



P=^exp[-ln- — =cosh (1 + 



x{(f)(q- < q >) 2 + iP(p- <p >) 2 - ( X ' + i-)[2(q- < q >)(p- < p >) - ih]}], (29) 



where 



Z = ^- x > 2 , y! = x-i\. (so) 

The density operator (29) is in a Gaussian form, as was expected from the initial form of the 
Wigner distribution function. While the density operator is expressed in terms of operators q and 
p, the Wigner distribution is a function of real variables q and p. When time t goes to infinity, the 
density operator approaches to 

P(°°) = ^=^ ex Ph^^ ln |^=^ ( 31 ) 
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where a = a pp (oo)a qq (oo) — a 2 „(oo) and [8]: 

(jpp(oo) = — ^ — - - ^ ((muj) 2 uj 2 D qq + (2A(A - fi) + uj 2 )D pp - 2mu 2 (X - n)D pq ), (32) 



2(mu) 2 \(\ 2 + u 2 - p 2 ) 
1 



O-na OO 



((mcu) 2 (2A(A + //) + u; 2 )/^ + tu 2 D pp + 2mtu 2 (X + (x)D pq ), (33) 



-(A + fj,)(mwYD qq + (A - (i)D„ + 2m(X 2 - (i*)D„). (34) 



m > 2m\(\ 2 + u 2 -v 2 ) 
In the particular case (7) 

ft , turf . . faw ftu; , . . . 

CT <ra °o = 75 coth 7777^ cr pp (oo) = — — coth — — , a pq {oo) = 35 

2muj 2kl 2 2kl 

and the asymptotic state is a Gibbs state (6): 

fru) 1.1., moo 2 a2 

PG (oo) = 2smh^ex P [-— (— p + —q )]. (36) 



Because of the presence of the exponential form in the density operator, the construction of the 
logarithmic density is straightforward. In view of the relations (25-27), the expectation value of 
the logarithmic density becomes 




(37) 



By putting hu = V5 — h/2, we finally get the entropy in a closed form: 



S(t) = k[{u+ l)ln(z/+ 1) -vlnv]. 



(38) 



It is worth noting that the entropy depends only upon the variance of the Wigner distribution. 
When time t -> oo, the function v goes to s = uj(d 2 2 /\ 2 - |di| 2 /(A 2 + ft 2 )) 1 / 2 /2ft - 1/2 and the 
entropy relaxes to its equilibrium value S(oo) = k[(s + 1) ln(s + 1) — sin s]. It should also be noted 
that the expression (38) has the same form as the entropy of a system of harmonic oscillators 
in thermal equilibrium. In the later case v represents, of course, the average of the number 
operator ||29|| . While the formal expression (38) for the entropy has a well-known appearance, the 
form of the function v displays clearly a specific feature of the present entropy. We see that the 
time dependence of the entropy is represented by the damping factor exp(— 2 At) and also by the 
oscillating function sin 2 (ftt). The entropy relaxes to its equilibrium value S(oo). 

5 Concluding remarks 

Recently we assist to a revival of interest in quantum Brownian motion as a paradigm of quantum 
open systems. There are many motivations. The possibility of preparing systems in macroscopic 
quantum states led to the problems of dissipation in tunneling and of loss of quantum coher- 
ence (decoherence). These problems are intimately related to the issue of quantum-to-classical 
transition. All of them point the necessity of a better understanding of open quantum systems 
and all requires the extension of the model of quantum Brownian motion. The Lindblad the- 
ory provides a selfconsistent treatment of damping as a possible extension of quantum mechanics 
to open systems. In the present paper we have studied the one- dimensional harmonic oscillator 
with dissipation within the framework of this theory. From the master equation of the damped 
quantum oscillator we have derived the corresponding Fokker-Planck equation in the Wigner W 
representation. The obtained equation describes an Ornstein-Uhlenbeck process. By using the 
Wang-Uhlenbeck method we have solved this equation for the Wigner function, subject to either 
the Gaussian type or the ^-function type of initial conditions and showed that the Wigner func- 
tions are two-dimensional Gaussians with different widths. Then we have obtained the density 
operator. The density operator in a Gaussian form is a function of q,p in addition to several 
time dependent factors. The explicit form of the density operator has been subsequently used to 
calculate the entropy. It relaxes to its equilibrium value. 
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